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The  translational  symmetry  breaking  of  a  crystal  at  its  surface  may  form  two-dimensional 
(2D)  electronic  states.  We  observed  one-dimensional  nonlinear  optical  edge  states  of  a  single 
atomic  membrane  of  molybdenum  disulfide  (MoS2),  a  transition  metal  dichalcogenide.  The 
electronic  structure  changes  at  the  edges  of  the  2D  crystal  result  in  strong  resonant  nonlinear 
optical  susceptibilities,  allowing  direct  optical  imaging  of  the  atomic  edges  and  boundaries  of  a 
2D  material.  Using  the  symmetry  of  the  nonlinear  optical  responses,  we  developed  a  nonlinear 
optical  imaging  technique  that  allows  rapid  and  all-optical  determination  of  the  crystal  orientations 
of  the  2D  material  at  a  large  scale.  Our  technique  provides  a  route  toward  understanding  and  making 
use  of  the  emerging  2D  materials  and  devices. 


The  structural  discontinuity  at  the  edges  and 
boundaries  of  2D  atomic  materials,  such 
as  graphene  and  transition  metal  dichal- 
cogenides,  leads  to  complex  interplay  between 
the  atomic  positions  and  the  electronic  structures. 
Subsequently,  the  atomic  edges  and  boundaries 
reconstruct  structurally  and  electronically  A  broad 
range  of  exceptional  physical  behaviors  and  ap¬ 
plications  including  widely  tunable  transport  gaps 
(7,  2),  unusual  magnetic  responses  (3-5),  and 
high-performance  nanoelectronics  (6,  7)  have  been 
discovered.  However,  experimental  observations 
of  these  ID  structures  have  been  limited  to  scan¬ 
ning  tunneling  microscopy  and  transmission  elec¬ 
tron  microscopy  Here,  we  studied  the  second-order 
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nonlinear  optics  on  the  ID  edges  and  boundaries 
of  hexagonal  molybdenum  disulfide  (MoS2) 
atomic  membranes.  The  broken  inversion  sym¬ 
metry  of  the  atomically  thin  monolayer  shows 
strong  second-harmonic  generation  (SHG),  in 
contrast  to  the  centrosymmetric  bulk  material, 
which  is  immune  to  the  second-order  nonlinear 
processes.  The  destructive  interference  and  anni¬ 
hilation  of  nonlinear  waves  from  neighboring 
atomic  membranes  reveals  the  few-atom-wide  line 
defects  that  stitch  different  crystal  grains  together, 
and  also  allows  the  mapping  of  crystal  grains  and 
grain  boundaries  over  large  areas.  Our  optical  im¬ 
aging  technique  enables  the  nonlinear  optical  de¬ 
tection  of  the  edge  states  at  the  atomic  edges  of  2D 
crystals  where  the  translational  symmetry  is  broken. 
The  observed  edge  resonance  of  SHG  clearly  in¬ 
dicates  the  electronic  structure  variation  at  the 
atomic  edges,  which  have  long  been  suspected  to 
be  the  active  sites  for  electrocatalytic  hydrogen 
evolution  (8). 

Unlike  gapless  graphene,  the  monolayer  form 
of  transition  metal  dichalcogenides  such  as  MoS2 
shows  a  direct  band  gap  at  visible  frequencies, 
making  them  emergent  semiconductors  for  nano¬ 


electronics  and  optoelectronics  involving  photo¬ 
voltaic  and/or  photoemission  processes  (9,  10). 
In  MoS2,  the  unique  local  orbital  properties  of  the 
heavy  transition  metal  atoms  and  broken  inver¬ 
sion  symmetry  of  the  monolayer  crystal  intro¬ 
duce  an  imbalanced  charge  carrier  distribution  in 
momentum  space,  giving  rise  to  a  novel  valley- 
specific  circular  dichroism  (11-12, 13).  Hexagonal 
bulk  MoS2  has  a  layered  structure  with  a  single 
layer  of  molybdenum  atoms  sandwiched  between 
two  layers  of  sulfur  atoms  in  a  trigonal  prismatic 
lattice.  Each  unit  cell  is  formed  by  two  mirrored 
sublattices,  making  the  bulk  crystal  centrosym¬ 
metric  and  prohibiting  the  second-order  nonlinear 
optical  processes  therein.  However,  in  a  monolayer 
consisting  of  one  layer  of  molybdenum  atoms  and 
two  sulfur  atom  layers,  or  few-layer  MoS2  with  odd- 
numbered  layers,  the  inversion  center  is  removed, 
resulting  in  a  strong  second-order  nonlinear  optical 
response.  The  origins  of  the  strong  nonlinear  opti¬ 
cal  responses  and  their  symmetry  properties  have 
been  discussed  (14-16).  The  observed  oscillatory 
nonlinearity  between  the  odd-numbered  and  even- 
numbered  layers  is  unique  to  all  layered  transi¬ 
tion  metal  dichalcogenides,  including  MoS2  and 
WSe2  (17). 

Because  SHG  is  strongly  affected  by  the  low¬ 
ered  symmetry  of  a  material,  the  structure  and 
symmetry  properties  of  MoS2  atomic  membranes 
can  be  probed  by  second-harmonic  emission,  which 
allows  us  to  demonstrate  the  nonlinear  optical  prop¬ 
erties  of  these  2D  materials.  Figure  1  shows  the 
linear  and  nonlinear  optical  images  of  a  continuous 
monolayer  MoS2  membrane  epitaxially  grown  by 
chemical  vapor  deposition  (CVD)  (18).  Oxide-on- 
silicon  substrates  with  oxide  films  285  nm  thick 
were  used  for  optimized  optical  contrast.  The  mono- 
layer  samples  show  high-quality  photoluminescence 
and  are  optically  uniform  over  large  areas  (Fig.  1  A). 
In  contrast,  the  SHG  image  (Fig.  IB)  reveals  the 
polycrystalline  nature  of  the  uniform  monolayer, 
and  the  average  grain  size  ranges  from  20  to  40  pm. 
Here,  the  incident  polarization  of  the  pump  beam 
is  along  the  y  axis  and  the  total  second-harmonic 
radiation  is  collected.  The  uniform  SHG  intensity 
within  each  grain  indicates  that  the  individual  grains 


Fig.  1.  All-optical  determination 
of  the  crystal  orientations  of  MoS2 
atomic  membranes.  (A)  Optical 
image  of  a  large  area  of  CVD-grown 
monolayer  MoS2.  (B)  SHG  image  of 
a  polycrystalline  monolayer  of  MoS2 
of  the  same  area.  The  grains  and 
grain  boundaries  are  clearly  revealed 
by  the  reduced  SHG  intensity  at  the 
boundaries.  The  crystals  are  connected 
by  faceted,  abrupt  grain  boundaries. 

The  scattered  bright  spots  are  from  the 
nanocrystals  produced  in  the  vapor  growth  process.  The  average  monolayer 
crystal  sizes  are  between  20  and  50  pm.  (C)  The  direct  crystal  orientation 
image  shows  the  crystal  orientations  of  the  irregular-shaped  polycrystalline 
aggregates.  (D)  Schematics  show  the  flakes  I  and  II  are  two  crystals  with 


opposite  orientations,  as  they  have  the  same  contrast  in  crystal  orientation 
maps  but  a  strongly  destructive  interference  boundary.  Crystals  II  and  III  show 
an  orientation  mismatch  of  -12°.  The  crystal  groups  a  and  b  show  two  cyclic 
twin  boundaries  with  60°  (a)  and  30°  (b),  respectively.  Scale  bar,  40  pm. 
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are  single  crystals.  At  the  grain  boundaries,  how¬ 
ever,  the  SHG  is  substantially  suppressed  as  a  re¬ 
sult  of  the  destructive  interference  and  annihilation 
of  the  nonlinear  waves  generated  from  the  neigh¬ 
boring  grains  with  different  orientations.  Although 
the  grain  boundaries  are  only  a  few  atom  sites  in 
width,  the  crystal  boundaries  are  clearly  seen. 
Moreover,  the  nonlinear  generation  not  only  re¬ 
veals  the  symmetry  properties  of  the  crystal,  but 
also  allows  the  rapid  mapping  of  crystal  orien¬ 
tations.  By  analyzing  the  polarized  components 
of  the  SHG  (see  supplementary  text  and  fig.  SI), 
we  can  map  the  crystal  orientation  of  the  poly¬ 
crystalline  atomic  membrane.  This  approach  al¬ 
lows  us  to  capture  a  complete  map  of  the  crystal 
grain  structures,  color-coded  according  to  crystal 
orientation  (Fig.  1C).  Large-area  CVD-grown  films 
clearly  have  complex  grain  geometries  and 
orientations. 

Because  of  the  threefold  rotation  symmetry 
of  monolayer  MoS2,  the  SHG  imaging  without 
phase  information  does  not  discriminate  between 
opposite  crystal  orientations.  Consequently,  it 


only  determines  grain  rotations  modulo  60°,  and 
the  measurable  difference  between  grain  orienta¬ 
tions  is  from  0°  to  30°.  For  example,  crystals  I 
and  II  have  no  contrast  in  the  color-coded  orien¬ 
tation  maps  (Fig.  1C),  but  the  clearly  visible  grain 
boundary  between  them  (Fig.  IB)  indicates  that 
the  two  crystal  grains  are  in  opposite  orientations. 
On  the  other  hand,  the  grain  boundary  between 
crystals  II  and  III  has  much  less  contrast  (Fig. 
IB),  but  the  two  crystals  are  distinctive  when 
looking  at  the  crystal  orientations  in  Fig.  1C 
(-12°  difference  in  orientation).  The  images  ob¬ 
tained  clearly  show  the  twin  grain  structures  and 
the  sixfold  symmetric  patterns  (Fig.  1,  B  and  C): 
Crystals  in  group  a  have  a  cyclic  twin  structure, 
whereas  each  crystal  in  group  b  instead  has  a 
-30°  difference  in  orientation  with  respect  to  its 
neighbors.  Such  a  characterization  process  con¬ 
ventionally  requires  diffraction-filtered  dark-field 
transmission  electron  microscopy  as  well  as  a  so¬ 
phisticated  and  time-consuming  analysis  proce¬ 
dure  (79).  Thus,  we  can  directly  visualize  the  grain 
structure,  providing  key  information  for  optimiz¬ 


Fig.  2.  SHG  image  of  discrete  triangular  islands  of  MoS2  crystals.  (A)  SHG  image  pumped  at 
1280  nm.  (B)  SHG  image  of  the  same  sample  pumped  at  1300  nm.  The  brim  of  crystal  shows  a  strong 
nonlinear  optical  (SHG)  edge  state.  Scale  bar,  20  pm.  (0  Cross  sections  [white  dashed  lines  in  (A)  and 
(B)]  compare  the  SHG  of  the  same  crystal  under  pump  wavelengths  of  1280  nm  (black)  and  1300  nm 
(red),  respectively.  The  enhanced  edge  responses  are  clearly  shown. 


Fig.  3.  Nonlinear  edge  resonance  of  MoS2  monolayer  membrane.  (A  to  F)  SHG  images  of  the 
same  monolayer  MoS2  sample  for  fundamental  wavelengths  of  1280  nm,  1290  nm,  1300  nm,  1310 
nm,  1320  nm,  and  1330  nm,  respectively.  Scale  bar,  20  pm.  (G)  Resonance  of  the  edge  nonlinear 
response  at  -0.947  eV  (1310  nm).  The  ratio  is  determined  by  the  integrated  SHG  with  the  areas 
circled  in  (A). 


ing  synthesis  strategies  for  large-area  CVD  growth. 
The  crystal  orientations  determined  by  our  non¬ 
linear  optical  method  correspond  well  to  trans¬ 
mission  electron  microscopy  observations  on  the 
same  sample.  Although  the  results  shown  here 
were  obtained  for  a  fundamental  wavelength  of 
1300  nm,  the  symmetry  and  the  polar  patterns  of 
SHG  are  ubiquitous  for  all  pump  wavelengths. 

Nonlinear  optical  susceptibilities  are  highly 
sensitive  to  extremely  small  changes  in  materials’ 
electronic  properties.  In  a  3D  crystal,  breaking 
the  translational  symmetry  at  an  interface  yields 
appreciable  electric  dipolar  contributions  to  the 
surface  optical  nonlinearity,  making  the  highly 
surface-  and  molecular-specific  second-harmonic 
spectroscopy  an  indispensable  tool  for  noninva- 
sive  study  of  surface  sciences  (20).  Not  only  does 
it  measure  dipolar  width  across  interfaces  (21),  it 
also  probes  the  real-time  dynamics  of  a  surface, 
such  as  atomic  reconstructions,  charge  transfer,  and 
molecular  conformational  transitions.  We  found 
that  the  broken  translational  symmetry  at  the  edge 
of  a  2D  crystal  potentially  creates  ID  atomic 
edges  with  localized  electronic  properties  that 
may  substantially  differ  from  the  2D  bulk.  Shown 
in  Fig.  2,  A  and  B,  are  SHG  images  of  discrete 
triangular  CVD-grown  monolayer  islands  at  pump 
wavelengths  of  1280  nm  and  1300  nm,  respec¬ 
tively.  Although  actually  rough  at  the  atomic  scale, 
the  edges  of  these  monolayer  islands  are  the  more 
energetically  favorable  termination  of  molybde¬ 
num  zigzag  edges,  and  they  appear  straight  and 
sharp  under  the  optical  microscope  (18).  The 
SHG  images  with  two  different  fundamental  wave¬ 
lengths  differ  at  these  atomic  edges.  When  pumped 
at  1300  nm,  a  brim  of  enhanced  SHG  is  observed 
at  the  edges  of  the  crystals  (Fig.  2B),  which  is  sub¬ 
stantially  stronger  than  that  of  the  central  region. 
At  the  comers,  the  SHG  is  even  stronger.  By 
comparison,  no  such  enhancement  is  observed  at 
the  edges  at  the  pump  wavelength  of  1280  nm. 
The  edge-enhanced  SHG  is  further  seen  in  the 
cross  sections  in  Fig.  2C.  Structure  and  field 
discontinuities  at  the  edges  of  these  2D  crystals 
can  potentially  lead  to  enhanced  in-plane  non¬ 
linear  polarizability,  like  the  surface-enhanced 
SHG  often  observed  on  metal  and  dielectric  sur¬ 
faces  (22,  23).  Such  enhancement,  however,  is 
equally  effective  for  different  pump  wavelengths 
as  long  as  there  are  no  additional  electronic  states 
located  at  the  edge.  The  stronger  second-harmonic 
edge  responses  pumped  at  1300  nm  versus  1280  nm 
are  therefore  the  results  of  the  electronic  structure 
variations  between  the  sharp  crystalline  edges 
and  the  bulk  monolayer  crystals.  Moreover,  as¬ 
suming  that  the  region  where  the  electronic  states 
are  affected  by  the  existence  of  the  atomic  edges 
is  about  the  same  as  the  thickness  (i.e.,  less  than 
1  nm),  it  requires  more  than  three  orders  of  mag¬ 
nitude  of  nonlinear  enhancement  to  exceed  the 
nonlinear  generation  in  the  central  region,  given 
an  illumination  pump  laser  spot  size  greater  than 
1  pm. 

To  further  probe  these  phenomena,  we  im¬ 
aged  the  SHG  at  the  edges  and  comers  of  the 
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Fig.  4.  Dispersion  of  the  nonlinear  per¬ 
mittivity  for  both  the  monolayer  MoS2 
and  the  nanoribbon  with  Mo-zigzag 
edges  derived  from  DFT.  The  monolayer 
(black)  shows  no  spectroscopic  features  (dis¬ 
persive  and  lossless)  at  the  low  photon  en¬ 
ergies,  and  the  features  between  1.5  and 
4.0  eV  are  associated  with  interference  be¬ 
tween  the  one-photon  and  two-photon 
processes.  The  nanoribbon  (red)  shows  a  two- 
photon  resonance  at  the  photon  energy 
~0.8  eV  due  to  the  subband  transitions 
from  the  valence  bands  to  the  isolated  edge 
states  originating  from  the  Mo-zigzag  edges. 

The  nanoribbon  considered  here  has  a  width 
of  5  unit  cells,  and  both  the  edges  are  termi¬ 
nated  by  Mo-zigzag  edges.  The  lattice  con¬ 
stants  were  assumed  to  be  the  same  as  the  monolayer,  and  the  Mo  edges  are  not  passivated  for 
simplicity.  The  overall  nonlinear  responses  are  averaged  over  the  entire  supercell  that  includes  both 
the  nanoribbon  and  the  vacuum;  the  enhancement  of  nonlinear  responses  could  therefore  be 
severely  underestimated. 


samples  at  different  fundamental  wavelengths. 
We  observed  a  resonance  wavelength  is  about 
1310  nm  (-0.95  eV),  a  slightly  higher  photon 
energy  than  the  half  transition  energy  of  the  A 
exciton  complex  (Fig.  3).  Quantum  mechanical 
expression  for  the  second-order  nonlinear  sus¬ 
ceptibility  of  a  material  involves  matrix  elements 
with  one  occupied  and  two  unoccupied  interme¬ 
diate  states.  When  there  are  additional  energy  levels 
nearly  coincident  with  virtual  transition  levels,  the 
coupling  between  the  radiation  and  the  nonlinear 
optical  susceptibility  becomes  large.  The  non¬ 
linear  optical  transition  matrix  elements  at  the 
edges  may  involve  additional  occupied  edge 
states  relative  to  the  crystalline  atomic  membrane, 
resulting  in  strongly  enhanced  second-order  non¬ 
linear  susceptibilities.  Density  functional  theory 
(DFT)  has  predicted  the  existence  of  these  edge 
states  and  provides  atomic-scale  insight  into  them 
(24);  such  ID  edges  have  been  identified  on  nano¬ 
particulate  MoS2  clusters  on  single-crystal  gold  by 
scanning  tunneling  microscopy  in  ultrahigh  vac¬ 
uum  (25).  The  existence  of  special  electronic  edge 
states  for  finite  low-dimensional  crystals  has  also 
been  reported  in  graphene  nanoribbons  (26,  27) 
and  at  the  end  of  ID  atomic  chains  (end  states) 
(28).  We  observed  the  strong  optical  nonlinear 
responses  of  such  edges  in  2D  MoS2  membranes, 
which  have  been  suspected  to  be  the  active  sites 
for  electrocatalytic  hydrogen  evolution  (8).  As  a 
control  experiment,  there  are  no  extraordinary 
photoluminescence  responses  observed  on  these 
edges  (fig.  S6).  It  is  also  worth  noting  that  there 
exist  prominent  excitons  in  thin-frlm  MoS2  (29) 
and  that  strongly  correlated  excitons,  especially 
those  with  high  binding  energies,  can  potentially 
have  an  impact  on  second-order  nonlinearities. 
Moreover,  the  local  doping  profile  at  the  bound¬ 
aries  and  edges  introduced  during  the  growth 
may  also  account  for  the  difference  at  the  edges 
and  the  center  of  the  monolayer  crystals. 

The  linear  and  nonlinear  optical  properties  of 
monolayer  MoS2  can  be  investigated  with  first- 


principles  calculations  (see  supplementary  text). 
The  probabilities  of  virtual  electronic  transitions 
can  be  derived  from  second-order  perturbation 
theory  (30).  As  shown  in  Fig.  4,  we  have  cal¬ 
culated  the  dispersion  of  the  nonlinear  optical 
susceptibilities  using  the  full  potential  linear  aug¬ 
mented  plane  wave  method  within  DFT.  Because 
MoS2  belongs  to  the  point  group  D3h,  there  is 
only  one  independent  component,  namely  = 
X$  =  X%1  =  xfl  =  X$-  In  Fig-  4,  the  tensor 
elements  % ^  of  monolayer  MoS2  and  of  a 
nanoribbon  with  zigzag  Mo  termination  are  com¬ 
pared.  The  monolayer  shows  no  spectroscopic 
features  below  the  band  gap  (-1.9  eV).  In  stark 
contrast,  the  zigzag  nanoribbon  shows  a  two- 
photon  resonance  at  -0.8  eV,  originating  from  the 
subband  transitions  from  the  valence  bands  to  the 
isolated  edge  states  of  the  Mo-zigzag  edges.  Such 
a  two-photon  nonlinear  optical  resonance  quali¬ 
tatively  explains  the  experimental  observations  at 
1300  nm.  In  our  calculations,  we  considered  an 
infinitely  long  nanoribbon  with  a  width  of  five 
unit  cells  with  Mo-zigzag  edge  terminations.  The 
lattice  constants  are  assumed  to  be  the  same  as 
the  monolayer,  and  the  Mo  edges  are  not  pas¬ 
sivated  for  simplicity.  The  overall  nonlinear  re¬ 
sponses  are  averaged  over  the  entire  supercell  that 
includes  both  the  vacuum  and  the  nanoribbon. 
Therefore,  the  enhancement  of  nonlinear  responses 
near  0.8  eV  could  be  severely  underestimated.  The 
spectroscopic  features  between  1  and  4  eV  are 
mainly  associated  with  interference  between  the 
one-photon  and  two-photon  processes,  which  are 
commonly  observed  in  nonlinear  semiconductors. 

Our  technique  has  revealed  second-order  non¬ 
linear  optics  at  the  boundaries  and  edges  of  the 
atomic  membranes  of  the  hexagonal  transition 
metal  dichalcogenide  MoS2,  with  the  strong  op¬ 
tical  resonance  indicating  a  distinct  edge  state  due 
to  the  translational  symmetry  breaking.  The  non¬ 
linear  response  reveals  the  structural  and  sym¬ 
metry  information  of  the  2D  atomic  monolayer 
based  on  interference  of  nonlinear  waves.  Such  a 


rapid  nonlinear  optical  imaging  technique  will 
allow  exploration  of  the  structural,  optical,  and 
electronic  properties  of  2D  atomic  layers  over  a 
large  area. 
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